An approximation for zero-balanced Appell function F\ near (1,1) 

D. Karp* 



Abstract. We suggest an approximation for the zero-balanced Appell hyper- 
geometric function F\ near the singular point (1, 1). Our approximation can 
be viewed as a generalization of Ramanujan's approximation for zero-balanced 
2F1 and is expressed in terms of 3^2- We find an error bound and prove some 
basic properties of the suggested approximation which reproduce the similar 
properties of the Appell function. Our approximation reduces to the approx- 
imation of Carlson-Gustafson when the Appell function reduces to the first 
incomplete elliptic integral. 



1. Introduction. The generalized hypergeometric function is defined by |XQ(, formula 4.1(1)] 

"j , \ _ y> {ai)k---{a P )k Z k 



pFq \h,...,b, 



00 



(bi) k ■ ■ ■ (b q ) k k\ ' 



where (a)o = 1, (a)k = a ( a + 1) ■ ■ ■ (a + k — 1), k = 1,2,..., is shifted factorial. This function is 
called zero-balanced if p = q + 1 and Y^h=\ a i = Yli=i h- 

Ramanujan (see O HI [5]) suggested the following approximations for zero-balanced 2 F\ an d 
3F2: 

B(a,b) 2 F 1 (a,b;a + b;x) = -ln(l-x)+j(a,b) + 0((l-x)ln(l-x)), x->l-, (2) 

where 

7(0,6) = 20(1) -#»)- 0(6), ^) = ^t4, (3) 

and 

n w£r 1 ^ ( Zm 03 |*) = -Mi-*)+£+ o(d - ,j m - ^ 1-, 

where 3? (03) > and 

L = 20(1) - tp{ax) - 0(a 2 ) + £ ^ ^ 7 ^ ■ 

fc(ai)fe(a 2 )fc 

These formulas have been generalized to q +iF q by N0rlund [T7], Saigo and Srivastava in |18j . 
Marichev and Kalla in [15] and Biihring in [7], see details in the survey paper by Biihring and 
Srivastava [8]. 

The Appell function F\ generalizes 2 F\ to two variables and is defined by [TO] : 



F 1 {a;f3 1 ,f3 2 ;'y;z 1 ,z 2 ) = 2^ r~\ rn ZlZ2 > ( 4 ) 



fe,n=0 
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for \zi\ < 1, | ^2 1 < 1 an d by analytic continuation for other values of z\, z%. An asymptotic 
expansion for F\ has been studied by Ferreira and Lopez in [11] in the neighborhood of infinity. 
This approximation can be converted into an approximation around (1,1) using the formula 

Fi(a; b, c; d; 1 — xz, 1 — yz) = x~ b y~ c F\ ( d — a; b, c;d;l , 1 ) . 

V xz v z ) 

It has been noticed by B.C. Carlson in [6] that the incomplete elliptic integral of the first kind is a 
particular case of F\ : 

F(X, k) = AFi(l/2;l/2,l/2;3/2;A 2 ,A; 2 A 2 ). (5) 

Carlson and Gustafson studied the asymptotic approximation for F(X,k) in [9|. Their expansion 
can be shown to be a particular case of the expansion for F\ given later in [11]. We will show below 
that both expansions (but not the error bounds!) can be obtained by simple rearrangement of ([3]) 
and use of known transformation formulas for F\. More precise approximations for F(X,k) which 
cannot be reduced to expansions from [11] have been given recently by S.M. Sitnik and the author 
in [H]. 

The purpose of this paper is to give an analogue of ([2]) for the "zero-balanced" Appell function 
Fi with 7 = a + ft + ft. Important properties of Fi are permutation symmetry 

F 1 {a;fa,fa;^;z 1 ,z 2 ) = Ft(a; fa, ft; 7; z 2 ,z x ), (6) 

reduction formulas 

F x (a;fa,fa;y,z, 1) = 2 F X {a, fa; 7; 1)2^1 (a, ft; 7 - fa\z), (7) 

F 1 (a;fa,fa;~/;z,z) = 2 F 1 (a,fa + ft;7;z), (8) 

and reduction formula ([5]). Our approximation reproduces the permutation symmetry ([6]), reduces 
to Ramanujan approximation given in ([2]) in cases given by ([7]) and (jHJ) and reproduces Carlson- 
Gustafson approximation for the values of parameters given in ([5]). 
Some new reduction formulas for F\ have been discovered in 1121 . 



2. Main results. To save space let us introduce the notation 

faM,b 2 ( x ^y) = B{a,bi + b 2 )Fi(a\bi,b 2 ;a + bi + b 2 ;x,y). (9) 
Our main approximation is given by 

' y - x 



( \ 1 1 1 / t , r x . h(y-x) f 1,1, 6 2 + l 

9aM,b2\x i y)=ai- +7(0.01 + 62) + 77— T \n \ sF2 9 h j_ h ju 1 

l — x (61 + o 2 )(l — x) \2,bi + b 2 + l 



1 — x 



(10) 



where 7(0, 61 + b 2 ) is defined in ([3]). The following theorem confirms that g a ,b!,b 2 is indeed a correct 
analogue of the righthand side of ([2]). 

Theorem 1 For < x < 1, < y < 1, a, b%, b 2 > 0: 

faMte(. x iV) = 9aMte( x 'y) + R aMte( x iV)i ( n ) 

with 

0<R aM}b2 (x,y) <r(l + a-aln(r)) = 0(rln(r)), (12) 

where in the last formula x,y — * 1, r = (1 — x)b\ + (1 — y)b 2 — > is the "rhombic" distance to 
x = y = 1 which is asymptotically equivalent to Euclidian distance. 

Corollary 1.1 Formulas (llip and (I12p imply in particular the inequality 

fa,b u b2 ( x , V) > 9a,b u b2 ( x , y) ( 13 ) 

for all x, y £ (0, 1). 
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Proof of Theorem [TJ A simple rearrangement of (j3|) gives 

a)fc(/V 

(rt) k k\ 



F 1 (a;P 1 ,p 2n ;z l ,z 2 )=f2 ^H^ 2F l (a + k,p 2n + k-z 2 )zt (14) 



k=0 

Suppose 7 = a + /3 2 , then 2 Fl in (fT4"j) is zero-balanced and we can apply [lOl formula 2.10(12)] 



r(r/ + /3) ^ (n!) 2 

(15) 



Fi(q; 0i, (3 2 ; a + (3 2 ; z 1 ,z 2 ) 
(n!) 2 ft! 



It gives 

r(a + ft) 

E ^ + k) M^ l)kZ ' [~ Ml ~ *0 + 2 ^(! + ") - ^0% + n) - ^(a + ft + w)](l - 3,)". (16) 

n,k=0 

Taking account of 

{a) k+n = (a) fc (Q + k) n = (a) n (a + n) k , 
the expression for Euler beta function 

r(a)r(ft) 

and the derivative formula 

2 Fi (a,b;c;x) = — 2 F 1 (a, b; c; x) = 2^ ry-^ ^{a) 2 F 1 (a, 6; c; ac), (17) 

fl fc=o ^ c ' fc 



identity (116p can be rewritten as: 

B(a, /3 2 )F 1 (a; /?i,/3 2 ; a + (3 2 ; z u z 2 ) = V {a ] n \^ )n (1 - z 2 ) n x 

x{[- ln(l - z 2 ) + 2-0(1 + n) - ip(a + n) - ip((3 2 + n)] 2 F 1 (a + n, ft; a; Z\) - 2 Fi(a + n, ft; a; z{)} . 

(18) 

Applying the transformation 

( 1 — w \ ( y — x 

Fi(a;h,b 2 ;a + bi + b 2 ;x,y) = ( — — J Fi f b\ + b 2 ; bx,a;a + bi + b 2 ; - - y 



to (|T6j) and (|T8|) in view of ([9]) gives 
/a,fei,fe 2 (x,y) = ( - — - ) x 



?l)jfe(6l + &2 + fc)n(a)n(- ln(l - y) + 2^(1 + n) - ip(a + n) - 0(&i + 62 + ft + n)) /y - x x '' 



ft!(n!) 2 (l-y)-« V*-* 



1 ~ y \ h j (a)n(fci + 



In — — + 2-^(1 + n) - V(a + n) - ^(61 + 62 + n) 

1 - y 



x 2 F X (h +b 2 + n, h; 61 + 6 2 ; |— |j (1 - y) n 
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Taking n = in the above formula and applying 



2-F1 h + 6 2 , h) h + 6 2 



y-x 

1 — X 



i-y 

1 — X 



-&1 



we get 



1 



1-1/ 
1 — x 



6i 



&1 + 62,61 
61 + 62 



y-x 

1 — X 



+ -R, (20) 



where it is clear from (1191) that 



R = 0((l-y)ln(l-y)), 



which is equivalent to the second formula in (|12|) . Formula (|20p can be easily put into a different 
form by differentiating the identity 







2^1(0, 


6; c; x) 


with respect to a: 






2 iV(a, 6; c; x) = - ln(l 


-x)(l- 


- x) c " a 


Hence: 








2 -Pi 


61 + 62,61 


y - x~ 


-a 


61 + 62 


1 — X 





\c— a— 6 



\c— a— b 77 



1 — X \ 1 — X I I — X x 

m ■ 



1-2// 1-2/ \i - y 



.Fx' 



"0,62 


y-x' 


61 + 62 


1 — X 



Since 



F (a, 6; c; z)| a=0 = 2^ irWk 



_y (PUk-l)\ j e _bz^ f 1,1,6 + 1 

1 V ;fc | a =0 fe=l V ^ 



2,c + 1 



(21) 



we will have 



iF,' 



"0,62 


y-x" 


61 + 62 


1 — X 



62(2/ - a;) 



i + 6 2 )(l-x) 3 2 V2,/5i + /3 2 + l 



l,l,/3 2 + l 



2/ ~ x 
1 — x 



In view of definition (j 1 Of) of Qafixfi^ix^y) formula (|20l) transforms into (j 1 1 1) . 

To estimate the remainder term we will use the ideas from p3] . An application of the integral 
representation |10l formula 5.8(5)] and a change of variable give (u = 1 — x, v = 1 — y): 



Fi (a; 61, 6 2 ; a + 61 + 6 2 ; 1 - u, 1 - v 

-*>l„,-&2 



u -»i v -02 Fl { &1 + &2 . b 62 . a + &1 + 62 . ! _ I ! _ I 



r(a + 6i + 6 2 ) /• t°- 1 (l + t)- a dt r(a + 6i + 6 2 ) 



r(a)r(6i + 6 2 ) 7 (l + nt) fo i(l + t;t)fe r(a)r(&i + 6 2 ) 



fa(t)h bl , b2 (u,v;t)dt. (22) 



where 



n-l 



fc=0 



/a(t) = 0(t a - 1 ), t^O /€^l,l- a> 



(23) 
(24) 



and 

h bub2 (u,v;t) 



(l + tu) b ^(l + tv) b 2 



El T\k-tkST~^ (&l)m(&2)fc-m m fc-m , ; , M , n 



k=0 



m=0 



m!(/c — m)! 
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h bub2 (u,v;t) =0(t~^- b2 ), t^oo => fteHo,^. 



(26) 



Spaces T and are defined in [14J. Basically, they mean nothing other than the asymptotic 
formulas satisfied by / and h, presented above. If a, b\, 62 are all positive conditions I and II from 
[14] are satisfied. 

Representation f|22j) is not precisely a Mellin convolution. However, if we approach the point 
u = v = (i.e. x = y = 1) along straight lines we can put u = 716, v = 72c, where 71 and 72 are 
positive constants and e — > 0. It this case 

h blyb2 {u, v > t) = hh ,62,71 ,72 

and (1221) takes the form of Mellin convolution. Since every point u, v lies on some straight line 
with endpoint (1, 1) and all our further speculations assume sufficiently small but fixed u, v there 
are always 71, 72 and e (of course non-unique) which are implied. Hence the theory from |14j can 
be applied. 
From 



OO , N 

k=n 



-l) n {a) n2 F 1 {a + n,l-l + n--l/t) (-!)"(«)„ f (1 - a 



t n+l n\ 



\n-l 



t n+1 (n - 1)1 J (1 + s/t) n+a 




ds 



it is obvious that sign(/ n ) = (— l) n . Similarly, from 

h blM , n (u,v;t) = f>l) fc i fc ( bl h {b2)k - m u m v k -™ 



k=n 



m=0 



ml(k — m)\ 



Y J (-l) k{ ^^2F l (b 1 , -k; l-b 2 -k; u/v). 



k=n 



it can be seen that sign(/i n ) = (—1)™. This shows that the remainder is always positive which 
implies in particular inequality (fT3|) . 

Now take n = 1 and apply |14} Theorem 4.3] which shows that the remainder has the form 
(since a = 0, b = 1 in terms of 



Ra,bxto{u,v) = J f a ,i(t)h bltb2t i(u,v;t)dt 




j.a—1 



(l + i) a t 



(l + ut) b i(l + vt) b z 



dt. (27) 



The bound for R a ,bi,b 2 ( u i v ) is based on the following lemma whose proof we postpone until the 
end of the proof of the theorem. 



Lemma 1 For all t G (0, 00) the inequalities 

- a/t 2 < f a>1 (t) < 0, 

- l/t < f ajl (t) < 0, 

- 1 < h blib2il (u,v;t) < 0, 
- t(ub\ + vb 2 ) < h bl fi 2l i(u, v; t) < 

hold true. 

The integral in (p7|) may be decomposed as follows 

1 l/r 00 

Ri = J fa,i(t)h bltb2t i(u,v;t)dt+ J fa,i(t)h buba ,i(u,v;t)dt+ J f a ,i(t)h blib2il (u,v)t)dt, 

1 l/r 



(28) 

(29) 
(30) 
(31) 
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where r can be any positive number (it is not needed that r < 1!). Set r = ub\ + vb 2 and use 
estimates ([29]) and ([3T]) in the first integral, ([28]) and ([3T]) in the second and ([28]) and ([30]) in the 
third. This gives the estimate (j!2p , □ 

Remark 1. We could use Proposition 3.1 from p3] to give an estimate for the error term. 
However, in our specific situation we are able to derive a much better bound based on Lemma Q] 
using the method of proof of this proposition but not it's statement. 

Proof of Lemma Q3 

(a) Inequality ([28]). Write / a ,i(i) = g a (t)/t 2 , where 

t a+1 

Then (|28p is equivalent to —a < g a (t) < 0. Clearly, g a (0) = 0. It is an easy exercise to check that 
g a (oo) = —a. If we prove that g' a (t) < we are done. Differentiating and multiplying both sides 
by (1 + t) a+1 we see that the required inequality takes the form 

(l + a)(l+t)t a <{l + t) a+1 + at a+1 & ^ +t)a+ . > 1 (l + x) a+1 > l + (l + o)x, 

t a (l + a + 1) 

where x = 1/t and the last inequality is the classical Bernoulli inequality valid for a > and x > — 1 
HQ formula 111(1.2)]. 

(b) Inequality (]29p is proved similarly but simpler. 

(c) Inequality ([30]) is obvious from the definition ([25]) of /i6 1) 6 2 (n, v; t). 

(d) To prove ([3T]) we again apply Bernoulli's inequality [16] formula 111(1.2)] in the form (pi, b 2 > 

0): 

(1 + tu)- bl > 1 - bitu, (1 + tu)- fc2 > 1 - b 2 tu. 
Multiplying these two inequalities we get the estimate 

1 - (i +tu) 4+^ < + v ' 2) - t2w6i&2 (32) 

which is even stronger than (|31|) . □ 

Remark 2. Application of (]32p instead of (|3ip in the proof of theorem [T] leads to an estimate 
of the remainder term R which is better than (j!2p . However, numerically it is only a very minor 
improvement, so we decided to keep the simpler estimate (]12p in the theorem. 

Remark 3. Representation ()10p also leads to the following observation: for general values 
of parameters there exists no approximation for f a ,b\to ^ n ^ e neighbourhood of (1,1) in terms of 
elementary functions. Indeed, let 

fa,b u b2( x >y) = H x ,y) 
as x,y — > 1 with an elementary h(x,y). Then from (jlip : 



h(x,y)+ln(l-x) = 2^(l)-V(a)-V(6i+6 2 ) + 7i-^T7^-T3i ? 2 ( ' ' ' ' ? " + ' 



7 — - ) +e(a:,y), 
1 — x 



(&i + &2)(l-a?) V 2 ^i + fe 2 + l 
and e(x,y) ^0 as x, y — > 1. Let x, ?/ — ► 1 along a straight line going through (1,1), so that 

(1 — y)/(l — x) = 7 = const. 

Then, due to 

(y - x)/(i - x) = 1 - 7, 

we have for the elementary hi(x, y) = h(x, y) + ln(l — x): 

hi(x,y) = d(i) +s(x,y), 
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where 

dfr) = w) - «.) - + ») + |^ s f 2 (^+ > + , 1 1 - 7 

Since y = 1 — (1 — x)7, we can write the above as 

/ii(x,7) = d(i) +e(x,i), 

where 7 G (0, 00) is arbitrary, but fixed. For x = 1 this gives ^1(1,7) = cZ('y) for all 7 £ (0, 00). 
Hence, a restriction of an elementary function h\ gives 3-F2 for all values of its argument in the 
range (—00, 1), which is impossible, and so h(x,y) cannot be an elementary function. 
Remark 4. Expansion formula (53)] can ce cast into the form 

r(bi + b 2 )T(a) / 71 72 

— — — -Fi a; 61, 6 2 ; a + 61 + 6 2 ; 1 , 1 

r(a + 61 + 62) V z z 



n-l 

£ 

fc=0 



IE + lQ g(^) IE 



+ -R n (a,&i,&2,7i>72;z), (33) 



Substituting x = 1 — 71/z, y = 1 — 72/z into (fT9|) we see that both (|33|) and (JT9J) are asymptotic 
expansions for \z\ — > 00 in the same asymptotic sequences z~ k log(z) and so their coefficients 
are the same. Hence, (|19j) can be viewed as a simpler form of formula (53)]. The appearance 
of the coefficients and Ej. is very different from that of the coefficients of (|19p and direct 
reduction is non-trivial. For instance, the first term of [111 formula (53)] reads (after some simple 
manipulations) (F = 2^1, M = (1 — y)/(l — x)): 



B(a,bi + b 2 )Fi(a;bi,b 2 ;a + h + b 2 ;x,y) = ip(l) - ip(a) 

1,62 + 1 



ln(l-^+ln(M)+^(l)-^(6 1 + 62) ( }fh>y 



+ 



61 + 62 
1 



61+62 + I 



1 - M 



\Mb 2 F' 



1,62 + 1 
61 + 62 + I 



1 - M 



+ b 1 F' 



1,62 

61 + 62 + I 



+ 61F 
1 - M 



1,62 

61+62 + I 



61 + 62 V 

Now using the relation |X0(, formula 2.8(36)] 

(c - a - b)F(a, b; c; z) - (c - a)F(a - 1, b; c; z) + 6(1 - z)F(a, 6 + 1; c; z) = 
we immediately get 



1 - M 
+ Ri. (34) 



(35) 



Mb 2 F 



1,62 + 1 
61 + 62 + I 



1 - M 



+ 61F 



1,62 

61 + 62 + I 



1 - M 



6i + 6 2 . 



Differentiating ([35]) with respect to a and putting a = we get: 

(c - 6 - 1)F'(1, 6; c; z) + 6(1 - z)F'(l, 6 + 1; c; z) = F(l, 6; c; 2) + (c - 1)F'(0, 6; c; z) - 1. 
Using (|2Tj) we see 



Mb 2 F' 

1,62 

61 + 62 + I 



1,62 + 1 
61 + 62 + I 



1-M 



+ 6iF' 



1,62 

61 + 62 + I 



1 - M 



1 - M 



62(61 + 6 2 )(1-M) 
(61 + 62 + I) 32 



1,1,62 + 1 
2, 61 + 6 2 + 2 



1 - M 



and 



S(a, 61 + b 2 )Fi(a; b±, 6 2 ; a + 61 + 6 2 ; s,?/) = In 
1 



1 



+ 



61 + 6 2 ' 



1,62 
61 + 62 + I 



1 - M 



+ 



1 — x 
6 2 (1 — M) 



+ 2^(1) - ip(a) - ipih + 6 2 ) 



(61 +62 + I) 



3F2 



1,1,62 + 1 
2, 61 + 6 2 + 2 



1-M 



1 



61 + 62 



+ Ri. 



(36) 
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Finally, (|36|) is reduced to (jlOp with the help of the following formula found at 
|http://functions.wolfram.com/07.27.17.0029.0l7] 

1 



iF 2 (a, b,c;a + 1, e; z) 



a — e + 1 



[a 2 -Fi(6, c;e;z) - (e - l) 3 F 2 (a, b, c; a + l,e- ljz)] 



Recalling that M = (1 — y)/(l — x) we get fjlOj) . The direct reduction for further terms is even 
more complicated. 

Theorem 2 The following properties are true: 

1. The function g is permutation symmetric: 

g a ,bi,b 2 {x,y) = 9a,b2,bi(y^ x )- ( 37 ) 

2. For y = 1 (and x = 1 due to (|37p ) the function g a ^ 6 2 (x, y) reduces to the Ramanujan's 
approximation: 



00,61,62 OM) = In -— — + 2^>(1) -V'(a) -ip(bi), 

g a , bl ,b 2 (l,y) = ln-— + 2tP(1) - ip(a) - i>(b 2 ). 
1-2/ 

5. For x = y the function g a , bl Mi x ->y) becomes the Ramanujan's approximation: 

g a ,b u b 2 {x,x) = In :j— !— + 2^(1) - ij>(a) - ip(bi + b 2 ). 

4- For the values of parameters a = b± = b 2 = 1/2 we have 

4 



(38a) 
(38b) 



01/2,1/2,1/2 



In- 



v / T 3 A2 + v / l- W 
which is the approximation of Carlson- GustaJ "son. 

Proof. To prove the first statement we need the following elementary lemma: 
Lemma 2 For b ^ 1 f/ie following relation holds true: 



(39) 



(40) 



3-^2 



1, b,c 

2, e 



z-l 



(l-z) b {c-e) 
c- 1 



"3-^2 



1, 6, e — c + 1 

2, e 



^ + 



( e -l)(l- z )(l- [1-zf- 1 ) 



For b = 1 it reduces to 

z 



3-^2 



1,1, c 

2,e 



z - 1 



( Z _l)( e _ c ) 
c- 1 



"3-^2 



1, l,e - c + 1 

2, e 



(c-l)(6-l)z 



: 1+ ^Ki^)ln 1 



(c-l)z 



1-z 



(41) 



(42) 



at 



http://functions.wolfram.com/07.27.03.0120.01): 



Proof. The proof is based on the following easily verifiable relation (which can be also found 

(43) 



3-^2 



1,6,C 

2,e 



e- 1 



2-^1 



6- l,c- 1 
e - 1 



(&-l)(c-l)z 

To prove (j41[) write this relation for zj[z — 1) in place of z, apply 



1 



2- r l 



6- l,c- 1 
e - 1 



z - 1 



(1 - z) b -\F x 



b — 1, e — c 
e- 1 



and substitute 2 -Fi from the right-hand side by 2 F\ expressed from (|43p . To prove (|42p let 6 tend 
to 1 and apply the L'Hopital rule. □ 



Combining ([32]) with the definition (fTU|) of g a ,b u b 2 { x iy) we immediately obtain ([37|) . 
Next we check the behavior of the function 9a,bi,b2(, x i 

y) on the sides of the square |x| < 1, 
\y\ < 1. Writing (|43p for z = 1 and using the Gauss formula for 2^1(1) we get 



3-^2 



l,6,c 
2,e 



e- 1 



(6-l)(c-l) 



2-^1 



b - l,c- 1 
e - 1 



v- 1 

e - 1 



(6-l)(c-l) 



r(e-l)r(e-fr-c+l) 
r(e - b)T(e - c) 



Now let 6 — > 1 and use the L'Hopital rule: 



3^2 



1,1, c 
2,e 

r(e) 



(e - l)r(e - 1) d r(e - b - c + 1) 



6=1 



(c - l)r(e - c) d& T(e - 6) 
-r(e - c)if>(e - c)r(e - 1) + T(e - l)^(e - l)r(e - c) _ (e - 1) 



~ ( C _i)r(e- C ) [r(e-i)P 

Substituting e = 61 + 62 + 1, c = 62 + 1 gives ([3"5j). 

Identity (f39|) is obvious from the definition (fTUj) of ga,bi,b 2 (x,y)- 
Finally, formula (I40p follows from the reduction formula 



(V(e-l)-^(e-c)). 



3-^2 



1,1,3/2 
2.2 



— In - + 



1 Vl^~z 



This completes the proof of the theorem. □ 
Corollary 2.1 For x, y — ► 1 



1 



+ 0(1). 



(44) 



1 — xy 

Proof. Assume first that a; and y approach (1, 1) in a way such (1 — y)/(l — x) stays bounded. 
We have 

1 



In 



1 — xy 



In 



1 



1 — x + x — xy 



In 



(l-x)(l + x±5f 



In 



1 



1 — x 



+ In 



1 



1 + x 



l-x 



Hence, 



In 



1 



1 — xy 



9a,bi,b 2 ( x ,y) 



In 



1 + x ±^y 

1 T X j_ 



7(0, 61 + 62 - 77 —7 V1 r3^2 

(01 + o 2 )(l - x) 



1,1,62 + 1 
2, 61 + 6 2 + 1 



1 — X 



O(l). 



If (1 — y)/(l — x) is unbounded, than exchange the roles of x and y and use f)3T|) . □ 

Finally, we remark that the authors of [U [2] consider monotonicity and ranges of the functions 

1 — 2-F1 (a, b; a + b; x) X2.F1 (a, b; a + 6; x) 



ln(l — x) 



ln(l/(l -x)) 



and 



B(a, 6)2-^1 (a, 6; a + b; x) + ln(l — x) 
for x 6 (0, 1). Our Corollary 12.11 shows that similar problems can be considered for the combinations 

1 - F 1 (a;Pi,f3 2 ;a + fa + fa;x,y) 



and 



for x, y G (0, 1). 



ln(l — xy) 



fa,/3i,/h(x,y) - hi 



1 



1 - xy 
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